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THE MASS FUNCTION OF ACTIVE BLACK HOLES IN THE LOCAL UNIVERSE 
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ABSTRACT 

We present the first measurement of the black hole (BH) mass function for broad-line active galaxies in the 
local Universe. Using the ~ 9000 broad-line active galaxies from the Fourth Data Release of the Sloan Digital Sky 
Survey, we construct a broad-line luminosity function that agrees very well with the local soft X-ray luminosity 
function. Using standard virial relations, we then convert observed broad-line luminosities and widths into BH 
masses. A mass function constructed in this way has the unique capability to probe the mass region < 10 6 M@, 
which, while insignificant in terms of total BH mass density, nevertheless may place important constraints on the 
mass distribution of seed BHs in the early Universe. The characteristic local active BH has a mass of - 1O 7 M 
radiating at 10% of the Eddington rate. The active fraction is a strong function of BH mass; at both higher and 
lower masses the active mass function falls more steeply than one would infer from the distribution of bulge 
luminosity. The deficit of local massive radiating BHs is a well-known phenomenon, while we present the first 
robust measurement of a decline in the space density of active BHs at low mass. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1 . THE LOCAL BLACK HOLE MASS FUNCTION 

There is strong evolution in both the number density and typ- 
ical luminosity of active galactic nuclei (AGNs) over cosmic 
time (e.g., Ueda et al. 2003; Richards et al. 2006), but from 
luminosity functions alone it is difficult to determine whether 
mass or luminosity evolution is the predominant agent of these 
changes. One thing we do know is that the growth of black 
holes (BHs) and galaxies are coordinated such that, in local 
spheroids, BH mass is strongly correlated with spheroid lumi- 
nosity (Marconi & Hunt 2003) and stellar velocity dispersion 
(the Mbh-c* relation; Gebhardt et al. 2000a; Ferrarese & Mer- 
ritt 2000; Tremaine et al. 2002; Barth et al. 2005). Charting the 
mass accretion history of the Universe thus may provide impor- 
tant insight into the growth of galaxies. It has become possible, 
using the Mbh-c* relation, to calibrate virial scaling relations 
between AGN luminosity and the size of the broad-line region 
(Gebhardt et al. 2000b; Ferrarese et al. 2001 ; Onken et al. 2004; 
Nelson et al. 2004; Greene & Ho 2006b); these techniques have 
been used to investigate BH mass functions at intermediate to 
high redshift (e.g., Vestergaard 2004; McLure & Dunlop 2004; 
Kollmeier et al. 2006), but not in a systematic way for nearby 
systems. A good measurement of the local active BH mass 
function provides an essential boundary condition for models 
of the evolution in active BH mass. Furthermore, we can probe 
significantly further down the mass function at the present day 
than at any other epoch. 

In fact, BH mass functions built from local broad-line AGN 
samples have the unique capability to probe BH masses Mbh< 
10 6 5 Mq. Mass functions derived from stellar velocity dis- 
persions (whether they be inactive [Yu & Tremaine 2002] or 
narrow-line active galaxies [Heckman et al. 2004]) are neces- 
sarily limited to the current spectral resolution limits of large- 
area spectroscopic surveys such as the Sloan Digital Sky Survey 
(SDSS; e.g., Bernardi et al. 2003). At the same time, the con- 



version from galaxy luminosity to BH mass is completely un- 
constrained at these low masses. Direct dynamical mass mea- 
surements in this mass range are beyond the spatial resolving 
power of current instrumentation for all but the nearest sys- 
tems. Thus broad-line AGNs currently provide the only means 
to systematically explore the BH mass function below 10 6 M© 
("intermediate-mass" BHs; e.g., Greene & Ho 2004). Although 
such objects constitute a negligible fraction of the present-day 
BH mass density, it is actually quite important to characterize 
the low-mass end of the local BH mass function. For one thing, 
it provides one of the only available observational constraints 
on models of the initial mass spectrum and halo occupation 
fraction of BH seeds in the early Universe (e.g., Volonteri et 
al. 2003). Furthermore, anisotropic gravitational radiation from 
unequal mass BH-BH mergers imparts a net linear angular mo- 
mentum or "kick" to the merger remnant with a velocity that 
may exceed the escape velocity of dwarf galaxies (e.g., Favata 
et al. 2004; Merritt et al. 2004). We may test this picture with 
observational constraints on the number of local dwarf galaxies 
that host BHs. 

Once we have characterized the zero-redshift broad-line BH 
mass function, we may investigate whether broad-line AGNs 
trace the same local population as samples selected by alter- 
nate means. For instance, we expect the soft X-ray luminosity 
to come from unobscured sources with broad lines, and thus 
we expect very similar luminosity functions for the two pop- 
ulations. Also, it would be instructive to compare our results 
to those of the complementary study of Heckman et al. (2004), 
which uses the Mbh - o~* relation to infer the BH mass func- 
tion for local narrow-line AGNs. Our methodology, based on 
AGN physics rather than indirectly on the Mbh _ c* relation, 
provides an important alternate measure of the local active BH 
mass density. Furthermore, we may compare the space density 
of narrow- and broad-line objects as a function of mass; a well- 
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matched comparison has never been performed in the literature 
before, and has interesting consequences for our understanding 
of AGN unification. Finally, as we shall see, the shape of the ac- 
tive and inactive (as inferred from galaxy luminosity functions) 
mass functions diverge, both at high and low mass. The mass- 
dependent active fraction contains useful information about the 
primary triggering mechanisms of active galaxies. By studying 
in detail the BHs that are radiating in the present-day Universe, 
we may hope to gain new insight into the strong evolution of 
BH growth with cosmic time. 

We briefly review the methodology behind BH mass esti- 
mates in §2, then present our sample selection and methodol- 
ogy in §3, derive the luminosity and mass functions in §4, and 
discuss the implications in §5. We summarize and conclude 
in §6. Throughout we assume the following cosmological pa- 
rameters to calculate distances: Ho = 100 h = 71 km s" 1 Mpc" 1 , 
fl m = 0.27, and n A = 0.75 (Spergel et al. 2003). 

2. BLACK HOLE MASSES FOR ACTIVE GALAXIES 

Compared to the dynamical masses used to calibrate the 
Mbh - c* relation, BH mass estimates using AGNs are rather 
crude, and some discussion of their validity is in order. So- 
called virial masses, Msu=fRv 2 /G, are estimated using a rela- 
tion between AGN luminosity and broad-line region radius (R; 
the radius-luminosity relation; e.g., Kaspi et al. 2005), com- 
bined with a measurement of the broad-line region (BLR) ve- 
locity dispersion (v). The factor / accounts for the unknown 
broad-line region geometry and here is assumed to be 0.75, cor- 
responding to a spherical broad-line region (Netzer 1990). 

Although relying on the physics of the BLR to estimate BH 
masses is fraught with danger (e.g., Krolik 2001), it is possi- 
ble to directly compare virial masses with the measurements 
of bulge stellar velocity dispersions, and calibrate them with 
the M B h - relation. Rather remarkably, those reverberation- 
mapped AGNs with er» measurements show good agreement 
between the two mass indicators (Gebhardt et al. 2000b; Fer- 
rarese et al. 2001; Onken et al. 2004; Nelson et al. 2004). Still, 
the total sample, across all above comparisons, consist of no 
more than 15 objects. For this reason Greene & Ho (2006b) di- 
rectly compared er* with virial masses for 88 AGNs. They find 
the virial masses provide a good estimate of BH with a scatter 
of 0.4 dex, but a number of potential systematic effects remain. 

The slope of the radius-luminosity relation is currently un- 
certain, primarily because it has been measured for only ~ 30 
objects spanning a limited luminosity range. Most recently, 
Bentz et al. (2006) find a shallower radius-luminosity slope 
(0.52 ± 0.04 rather than the 0.64 ± 0.02 assumed here), which 
would lead to a narrower distribution in BH mass. For the typ- 
ical Eddington ratio of our sources, Lboi/£Edd=0.1, the Bentz 
et al. formalism shifts the BH masses upwards by ~ 0.3 dex 
for a 10 5 Mq BH and downwards by ~ 0.2 dex for a 10 9 M Q 
BH. However, as noted by Netzer & Trakhtenbrot (2007), if 
all objects included in the Kaspi et al. sample are included in 
the Bentz et al. analysis, the final slope is consistent with the 
Kaspi et al. result. Larger reverberation-mapped samples are 
urgently needed, but in the meantime we continue to use the 
steeper slope. The other outstanding uncertainty is related to 
the (unknown) geometry of the BLR, which translates into a 
different pre-factor /. Both Onken et al. (2004) and Greene & 
Ho (2006b) find that virial masses calculated assuming a spher- 
ical BLR (as we have done here) are too low by a factor of 
~ 1 .6, corresponding to a net increase in all BH masses of 0.2 
dex. However, Collin et al. (2006) present some evidence that 



the BLR geometry is luminosity or Eddington-ratio dependent. 
As yet, there are not sufficient statistics to define a luminosity- 
dependent / factor, and so we prefer to apply no correction at 
this time. 

3. SAMPLE SELECTION 

The sample selection follows the technique of Greene & Ho 
(2004), which we summarize briefly here. Our parent sample 
is drawn from the Fourth Data Release (DR4) of the SDSS 
(York et al. 2002; Adelman-McCarthy et al. 2006), which in- 
cludes spectroscopy of 565,715 objects classified as galaxies 
and 67,382 objects classified as AGNs with redshifts z < 2.1. 
We analyze all spectra with z < 0.352, although we remove 
those galaxies for which > 20% of the pixels between 6400- 
6700 A have been flagged as bad by the SDSS pipeline; this 
results in a total of 544,127 spectra (including multiple obser- 
vations of the same object). Following Greene & Ho (2004; 
see also Ho et al. 1997a), we begin by modeling and removing 
the stellar continuum, which may mask or mimic the presence 
of broad emission lines that we use as the signature of an ac- 
creting BH. We use the principal component analysis (PCA) 
method of Hao et al. (2005a) that encapsulates the variance 
in SDSS absorption-line galaxies in a set of eight orthogonal 
eigenspectra. The galaxy continua are modeled as linear com- 
binations of these eigenspectra, a power-law to represent the 
AGN continuum, and a possible A-star component. We have 
not attempted to model the broad Fe II emission that is present 
throughout most of the optical spectrum (e.g., Boroson & Green 
1992; Greene & Ho 2006a); while this is important for mea- 
surements of H/3 and [O III], the impact on the Ha region is 
negligible. 

After galaxy subtraction, in order to limit the number of po- 
tential candidates, we employ an algorithm that identifies po- 
tential broad emission in the Ha region, based on an excess 
variance in the Ha region compared with the featureless con- 
tinuum. The following analysis is performed over the spectral 
region 6400-6700 A. Narrow emission-line flux is removed in 
an iterative scheme, in which points with fluxes > 4 a (as mea- 
sured in the entire spectral region) are replaced with the median 
value in the region, until there is < 1 % decrease in <r, in prac- 
tice taking no more than 5 iterations. This method of narrow 
emission line removal is conservative, in the sense that some 
flux from very strong narrow emission lines will remain. How- 
ever, without a full spectral decomposition of the region, which 
we perform below, a single-Gaussian fit to the narrow emission 
lines may easily remove a relatively narrow broad component. 
We then smooth the spectrum with a 10-pixel boxcar to remove 
high-frequency noise and increase our sensitivity to real broad 
components. Those sources with a root-mean-square (rms) de- 
viation that is 50% higher in the 6400-6700 A region than in the 
regions 5500-6280A (excluding 20 A around the NaD A5895 
doublet) and 6900-7300 A are kept, for a total of 31,750 spec- 
tra. 

With this smaller subset of sources we perform more detailed 
modeling of the Ha+[N II] AA6548, 6583 region. Careful mod- 
eling of the narrow lines is crucial, since velocity structure in 
the narrow-line region can mimic broad lines (see additional 
arguments in Ho et al. 1997c and details of the methodology 
in Greene & Ho 2004, 2005a,b). Whenever possible, we build 
a multi-component Gaussian model for the narrow lines using 
the nearby [S II] AA6716, 6731 lines that empirically are found 
to provide a good model for the narrow Ha+[N II] (Filippenko 
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FIG. 1. — Distributions of broad Ln a (left) and broad FWHMh q (right) for the entire sample of broad-line AGNs. The shaded histograms indicate the distribution 
of those objects targeted within the SDSS main galaxy sample. Note the dominance of objects with FWHMh q < 2000 km s _1 , the canonical division between 
classical and "narrow-line" Seyfert 1 galaxies. 



& Sargent 1988; Ho et al. 1997c), and in cases without [S II] 
emission, we resort to using a model based on the core of the 
[O III] A5007 line. The relative wavelengths of Ha and the 
[N II] doublet are fixed to laboratory values, as are the relative 
strength of the [N II] lines, and only the relative amplitudes of 
the Ha and [N II] lines are allowed to vary. We adopt the em- 
pirical criterion of Hao et al. (2005a) to determine whether an 
additional component is required: any component that results 
in a 20% decrease in \ 2 is deemed statistically justified. Broad 
lines are often lumpy and asymmetric, so we model the broad 
emission with multiple Gaussians, but attach no physical sig- 
nificance to any given component. All objects with broad Ha 
components apart from those in the narrow-line model are re- 
tained. 

Table 1. The Sample 



Name 


z 


FWHM Hq 


log L Hct 


log M BH 


(i) 


(2) 


(3) 


(4) 


(5) 


J000000. 47-002703. 9 


0.250 


2010 


41.65 


6.7 


J000048. 15-095404.0 


0.206 


2470 


41.61 


6.9 


J000102. 19-102326. 8 


0.294 


5330 


42.75 


8.2 



Note. — Col. (1): Name. Col. (2): Redshift. Col. (3): FWHM Ho 
(kms -1 ). Col. (4): L Ha (ergs s" 1 ). Col. (5): M BH = (2-0 X 10 6 

M (Lhc/10 42 ergs s" 1 )° ^±0.02 (FWH M Hq /10 3 km s -l)2-06±0.06 
(Greene & Ho 2005b). Table 1 is available in its entirety via the link to 
the machine-readable table above. The following is included only as a 
guide to content and presentation. 

Our goal is to limit our sample to those targets for whom 
we may estimate a reliable BH mass from broad Ha. We 
thus very aggressively flag all suspicious objects based on their 
line strength, by applying a combined cut on total Ha flux 
normalized to the rms deviations in the continuum-subtracted 



spectrum, and the Ha equivalent width (EW). We require that 
/ Ha /rms > 200 and EW(Ha) > 15 A. The EW measurement 
is based on continuum measured from the best-fit PCA galaxy 
model in the spectral regions 6570-6590 A and 6480-6540 A. 
In detail, our thresholds are selected based on simulations, pre- 
sented in Appendix A, which demonstrate that for the S/N ra- 
tios typical in the SDSS, BH mass estimates are highly uncer- 
tain (> 1 dex) and significantly biased for lower Ha fluxes and 
EWs. Hereafter, this combination of cuts will be referred to as 
the "detection threshold." Following the automated flagging, 
some objects at the low-mass end are still questionable when 
examined by eye. They are flagged as well. A total of 11,428 
objects are removed at this stage. For comparison with previous 
work, we also build a more inclusive luminosity function that 
does not reject the flagged objects. 

There are a number of scientifically interesting objects that 
fall below the detection threshold. There are strongly star- 
forming galaxies with a broad Ha base (see also Hao et al. 
2005a) that are probably Wolf-Rayet galaxies (Ho et al. 1997c). 
In other cases, we are probably seeing a small, scattered com- 
ponent of broad emission in what is predominantly a narrow- 
line AGN. In this case, the Ha luminosity does not correlate 
with the broad-line region size. Finally, if the BH is radiat- 
ing at a very low fraction of its Eddington luminosity then it is 
thought that the standard optically thick, geometrically thin ac- 
cretion disk (Shakura & Sunyaev 1973) is replaced by a radia- 
tively inefficient accretion mode (e.g., Quataert 2001; Narayan 
2006 and references therein). Although no such objects are in- 
cluded in the reverberation-mapped samples used to calibrate 
the radius-luminosity relation (e.g., Kaspi et al. 2005) it is likely 
that the broad-line region structure is quite different in this 
regime. While we do not have sufficient S/N in single objects 
to accurately measure BH masses, it is quite possible that in 
stacked spectra we would be able to uncover very low-contrast 
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FIG. 2. — Distributions of BH mass (left) and Lbol/^-Edd ("g^O for the entire sample of broad-line AGNs. Mbh is calculated from Ln a and FWHMho 
as shown above, using the formalism of Greene & Ho (2005b), and the Eddington ratio is derived assuming an average bolometric correction of L bol = 
2.34 X 10 44 (Z,Ha/10 42 ) ' 86 ergs s -1 (see text). As in the previous figure, the galaxies targeted in the SDSS main sample are indicated with the shaded histogram. 
Note that while the galaxy-selected AGNs are uniformly at low luminosities, they span the entire range of Mbh in the sample. Massive BHs are predominantly in 
low-accretion states at the present time. 



broad lines. We defer such analysis to a later work. 

Finally, we note that narrow-line AGN samples are typically 
selected using their location in diagnostic diagrams (Baldwin et 
al. 1981; Veilleux & Osterbrock 1987; Ho et al. 1997a; Kauff- 
mann et al. 2003; Hao et al. 2005a) that discriminate the shape 
of the ionizing continuum by the relative strengths of various 
prominent narrow emission lines. Our selection technique re- 
lies on the presence of broad emission lines, and so is unbiased 
with respect to line ratios. This is particularly useful for select- 
ing relatively metal-poor AGN hosts, since the typically high 
ratio of [N II]/Ha may be significantly reduced. This has been 
observed in NGC 4395 (e.g., Kraemer et al. 1999), with low 
significance in the original Greene & Ho (2004) sample, and 
quite strikingly in our new sample of Type 2 low-mass BHs 
(A. J. Barth et al., in preparation; see also Groves et al. 2006). 
In contrast to standard narrow-line selections, therefore, our se- 
lection is relatively unbiased against AGNs in low-metallicity 
(and thus typically low-mass; e.g., Tremonti et al. 2004) galax- 
ies in this manner. 

Our final sample (with duplicates removed) of broad-line 
AGNs is composed of 8728 objects. The Princeton spectral re- 
ductions include single-Gaussian fits to narrow emission lines, 
and flags for all objects with an additional broad component. 
Approximately 30% of the broad-line objects in their sample 
are rejected by our algorithms, some (35%) with the initial 
Ha-finding algorithm, and the rest either because of a lack of 
broad Ha in the multi-Gaussian fit, or because they fall below 
our detection threshold. We have examined a large number of 
these objects visually, and we feel justified in rejecting them. 
One contaminant in particular consists of objects with multi- 
component narrow lines, which confuse algorithms that fit only 
single Gaussians to the narrow lines. We also compare with 
the Hao et al. (2005a) list of broad-line AGNs. Only 53% of 



their broad-line sample is included in our final sample; 17% 
are rejected by our initial broad-line search algorithm, 1 1% are 
deemed to have no broad Ha in more detailed fitting, and 19% 
are rejected by our detection threshold. 

We have measured the Ha luminosities (with no internal red- 
dening correction applied) and FWHM Ha for our entire sam- 
ple of broad-line AGNs (Fig. 1). The FWHM Hq are measured 
from the multi-component Gaussian models (as in Greene & Ho 
2004). We can use these measured quantities and the formal- 
ism presented in Greene & Ho (2005b) to derive BH masses 
for the sample. We also use the measured Ha luminosities 
and the inferred BH masses to estimate Eddington ratios for 
the sample, L bol /L EAd , where L E dd = 1.26 x 10 38 (M BH /M Q ) 
ergs s" 1 . Assuming that Lbol = 9.8 L5100 (McLure & Dun- 
lop 2004), in terms of L Ha our bolometric correction is Lb i = 
2.34 x 10 44 (L H a/10 42 ) ' 86 ergs s" 1 (Greene & Ho 2005b). We 
adopt this bolometric correction throughout. The correspond- 
ing distribution of BH masses and Eddington ratios is shown in 
Figure 2. At face value the mass function is strongly peaked 
at 10 7 Mq, while the typical Eddington ratio of our sample is 
about a tenth of the Eddington limit. However, there are severe 
selection effects that cause us to lose significant numbers of 
sources, and these depend on mass, Eddington ratio, and red- 
shift. The strong importance of selection effects is shown in 
Figure 3, where the distribution of BH mass in the sample is 
shown as a function of both redshift and Eddington ratio. The 
most striking trends in this figure are the decrease of character- 
istic BH mass as one moves to higher Eddington ratio at a given 
redshift, and the increase of characteristic mass as one moves to 
higher redshift at a given Eddington ratio. Both of these trends 
are driven by the magnitude limits of the SDSS. 
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FIG. 3. — Distribution of Mbh with redshift is shown in small points for the 
entire sample. In large symbols we represent the mean mass in three different 
Eddington ratio bins, with 1 a contors bracketing each bin. The blue crosses 
and dash-dotted lines show objects with Lbol/^Edd< 0.01, the black triangles 
and dashed lines show 0.01 <£bol/^-Edd< 0.1, and the red circles and solid 
lines are objects with A,ol/^Edd> 0.1. 
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FIG. 4. — Volume-weighted broad Ha luminosity function in bins of 0.25 
dex (# Mpc~ 3 log Lh q ~'). The maximum volume is calculated based on both 
the photometric and spectroscopic limits of the survey and our search algo- 
rithm (see text). The error-bars represent the Poisson errors in each bin. The 
solid line represents our maximally inclusive sample, in which objects with 
less reliable broad Ha are included. Our best-fit double power law is shown as 
the dashed line. The inset panels show the luminosity functions for objects tar- 
geted as galaxies or AGNs, respectively, primarily based on a color selection. 



4. LUMINOSITY AND MASS FUNCTIONS 
4. 1 . Broad Ha Luminosity Function 



We begin by calculating a traditional luminosity function 
&(Lho) for our broad-line active galaxy sample. While the 
luminosity function is difficult to interpret directly, in so far 
as it is a convolution of the BH mass and Eddington ratio 
distributions, nevertheless it has the benefit of being a di- 
rectly observable quantity. BH masses, on the other hand, 
are potentially subject to considerable systematic uncertain- 
ties (e.g., Greene & Ho 2006b; Collin et al. 2006). Further- 
more, the luminosity function is readily compared with var- 
ious probes of the active galaxy population, such as the X- 
ray luminosity function for AGNs. Following common prac- 
tice, we present the differential luminosity function, the num- 
ber of AGNs per unit volume per logarithmic luminosity inter- 
val, $(Lea) = (LHa/log 10 e)$(Lii a ). Since galaxies of different 
magnitude are observable to varying distance within the sam- 
ple, a proper representation of the distribution of AGN luminos- 
ity (or mass) requires that we somehow account for our variable 
sensitivity prior to comparing their distributions directly. We 
primarily use the classical V /V mlix weighting method (Schmidt 
1968; Huchra & Sargent 1973; Condon 1989; Ulvestad & Ho 
2001), but in Appendix B we present a non-parametric maxi- 
mum likelihood luminosity function as well (e.g., Efstathiou et 
al. 1988). 
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FIG. 5. — Broad Ha luminosity function calculated using the maximum 
likelihood formalism (solid line; see Appendix B for details). The V/Vmax lu- 
minosity function is shown in solid symbols as in Figure 4. Also shown are 
the single (dashed line) and double (dot-dashed line) power-law fits to the Hao 
et al. (2005b) maximum likelihood total Ha luminosity function. 



In the V /Vmax formalism, our sensitivity is characterized as 
the maximum volume to which each source would be included 
in our final sample. In calculating V max , we must account for 
both the magnitude limit of the SDSS and the S/N dependence 
of our Ha detection procedure. The appropriate magnitude 
limit will depend on whether a given object was targeted as a 
galaxy or quasi-stellar object (QSO). The main galaxy sample 
comprises spatially resolved (i.e., non-stellar) targets above a 
limiting r-band Petrosian magnitude of 17.77 mag (Strauss et 
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al. 2002). QSO targets, in contrast, are selected with a rather 
complicated set of color-selection criteria within a magnitude 
range of 15.0 < i < 19.1 (Richards et al. 2002; Stoughton et 
al. 2002). The maximum flux limit mitigates bleeding be- 
tween fibers. A different color criterion is used to select high- 
redshift QSO candidates, and these have a limiting magnitude 
of i = 20.2; a small number (237) of our final targets fall into this 
category. We note that it is possible for objects to have been tar- 
geted as either a QSO or a galaxy; these objects are considered 
part of the galaxy sample. Now, in addition to these two pri- 
mary surveys, some "serendipitous" sources are targeted, due 
to a detection either in the ROSAT All-Sky Survey (Voges et al. 
1999), or by the FIRST radio survey (Becker et al. 1995; see 
also Richards et al. 2002). We have 154 such objects; they do 
not constitute a complete sample, nor do they have a uniform 
detection limit. We therefore assign them a maximum volume 
corresponding to their observed volume (i.e., take V max = V), 
and then verify that our conclusions are unchanged when we 
exclude them from the sample. 
Table 2. Broad Ha Luminosity Function 







log (L Ho ) 


log 


(1) 


(2) 



39.75 


-5 


91 + 0.24 
M -0.60 


40.00 


-5 


g5 +0.14 


40.25 


-5 


, q +0.10 
° -0.13 


40.50 


-5 


ln +o.08 

-0.10 


40.75 


-4 


95+0.07 
a -0.09 


41.00 


-5 


06 +0.04 
U -0.04 


41.25 


-5 


14 +0.03 
i4 -0.03 


41.50 


-5 


„ o +0.02 
z -0.03 


41.75 


-5 


42 +0.02 
* -0.02 


42.00 


-5 


6 „ + 0.03 
OD -0.03 


42.25 


-5 


„. + 0.03 
ao -0.03 


42.50 


-6 


9B+0.03 


42.75 


-6 


.,, + 0.04 
' J -0.04 


43.00 


-7 


„„+0.04 
^ -0.04 


43.25 


-7 


= 8 +0.07 



trast of the broad line, the galaxy subtraction, the strength of 
the narrow lines, etc. In principle, it is also possible that the 
galaxy dilution increases as the physical radius subtended by 
the fixed SDSS fiber increases, although it will be mitigated to 
some degree by surface brightness dimming, and we have ne- 
glected this effect. To model our incompleteness, we generate 
artificial spectra over a grid of redshifts (S/N) between the ob- 
served redshift and either the photometric redshift limit or the 
detectability threshold redshift, whichever is smaller. The S/N 
at each interval is calculated based on the fiber magnitude of 
the source at the new redshift and an empirical relation between 
fiber magnitude and S/N we have derived from the SDSS data. 
If the interval is smaller than Sz = 0.05, we simply bisect the 
redshift interval, whereas for larger redshift ranges we generate 
artificial spectra at z max _ 0.05 and the center of the redshift in- 
terval. We generate three spectra at each redshift, with a range 
of S/N ±20% of the calculated value. When the calculated S/N 
is < 5, the limiting redshift is automatically reached, as such 
spectra are simply unusable. We subject the remaining spectra 
to our full selection algorithm. As long as two of the three ar- 
tificial spectra are recognized as broad-line AGNs, we continue 
to the next redshift bin. In practice, only ~ 20% of our V max 
calculations are lowered by these simulations, suggesting that 
our detection threshold is reasonable. 

Using the V /V max weights derived above, the luminosity 
function ($) and uncertainty per bin (a) are simply calculated: 



'max / [ 



1/2 



(1) 



Note. — Col. (1): Log 
L Ha (ergs s- 1 ). Co^ (2): 
Space density Mpc~ ^ log 
Lji a —1 ) as a function of 
broad Ha luminosity. 

Based on the target selection of each object, it is straight- 
forward to calculate the maximum volume to which the SDSS 
would have spectroscopically targeted a given object, with ap- 
proximate k-corrections based on the observed colors. In de- 
tail, of course, the color (and thus target selection) of a source 
is redshift-dependent, as the spatially resolved galaxy light will 
dim more rapidly than the unresolved AGN light. This is diffi- 
cult to model without a good separation of host and AGN light. 
Nevertheless, since the majority of the galaxy targets are very 
close to their limiting magnitude, this should not affect the re- 
sults in practice. 

With increasing redshift and commensurately degraded S/N, 
it becomes more and more challenging to detect broad Ha. At 
a certain point, the object reaches our imposed detectability 
threshold, and is removed from the sample. In detail, how- 
ever, our algorithm may lose the ability to detect the object 
even before this limit is reached. This depends on the con- 



they are shown in Table 2 and Figure 4. The survey volume is 
determined assuming an angular coverage of 4783 deg 2 of the 
DR4 spectroscopic survejQ. We compute the number of AGNs 
per unit volume, per unit logarithmic luminosity [$(Lhq)] in 
bins of 0.25 dex in Ha luminosity, from 10 40 - 10 43 25 ergs s . 
In the insets are plotted the luminosity functions for objects tar- 
geted as galaxies (left) and QSOs (right). The "galaxies" are 
strongly peaked at the luminosity of a 10 7 M Q BH radiating 
at 10% of Eddington. At higher luminosity, the QSOs tend to 
dominate the light, while at the lowest luminosities, there are 
equal numbers of each type, depending on whether the source is 
a low-mass BH at higher Eddington ratio (QSO) or a more mas- 
sive and less active source. Note that the decline at luminosities 
below < 10 40 75 ergs s" 1 is due to incompleteness. For instance, 
the Palomar spectroscopic survey of nearby galaxies, with sig- 
nificantly higher sensitivity to low-level broad emission, finds a 
luminosity function that continues to rise to significantly lower 
luminosities (Ulvestad & Ho 2001; Ho 2004). For this reason, 
we also calculate a maximally inclusive luminosity function, 
retaining all nagged objects. That is shown as a solid line in 
Figure 4. 



^ |http : / /www, sdss . org/dr4/| 
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FIG. 6. — Conversion of our observed luminosity function to an X-ray lumi- 
nosity function using the Strateva et al. (2005) relation between UV luminosity 
and a ox . Solid line is the soft X-ray, unabsorbed luminosity function presented 
in Hasinger et al. (2005). The points represent the median luminosity function 
of 1000 Monte Carlo realizations of the optical to X-ray conversions, which 
account for its significant scatter. The dotted lines bracket 68% of the resulting 
distribution in a given luminosity bin. 

For comparison with other samples it is useful to parameter- 
ize the observed luminosity function. Following previous work, 
we consider a Schechter (1976) function and a double power 
law (e.g., Croom et al. 2004; Hao et al. 2005b), 

P(L*)/L. 



<j>{L) = 



(2) 



(L/L*r + (L/L^ , 

where </>*(L*), L„, a and (3 are free parameters. We use a x 2 
minimization fitting procedure that accounts for the asymmetric 
errors in space density, and we fit all objects in our final sam- 
ple with luminosities > 10 40 ergs s" 1 (for this luminosity range 
the results do not change if we fit the inclusive sample). The 
Schechter fit has a reduced \ 1= H> while the double power- 
law fit has a reduced xl= 1-3. Therefore, we prefer the lat- 
ter, which yields values of </>(L*) = (2.0 ±0.4) x 10" 6 Mpc" 3 , 
L* = 10 420±01 ergs s" 1 , a = 1.29 ±0.09, and (3 = 2.82 ± 0.07. 
This fit is shown as a dashed line in Figure 4. 



4.2. Comparison with Hao et al. (2005b) 

Hao et al. (2005b) used a smaller subsample of the SDSS 
to compute emission-line luminosity functions for broad- and 
narrow-line AGNs. They used very comparable selection crite- 
ria to select their broad-line subsample, and thus it is important 
to compare our results with theirs. In order to compare directly 
with their work we construct an additional luminosity function 
using the maximum likelihood formalism (for details, see Ap- 
pendix B). At high luminosity, Hao et al. have limited statistics 
(as demonstrated by the divergence of their single and double 
power-law fits in Fig. 5), but the double power-law fit does pro- 
vide a better fit to their highest luminosity bins, consistent with 
our result. At lower luminosity, our function diverges from that 
of Hao et al. We believe this discrepancy may be attributed to 
three factors. 

Firstly, we consider flux only from broad Ha, while Hao et 
al. include both narrow and broad emission. This has no effect 



at high luminosity, where the narrow-line luminosity constitutes 
a negligible fraction of the total, but at lower luminosities, ob- 
jects with high-EW narrow lines will be preferentially shifted 
into higher luminosity bins in the Hao et al. function. Secondly, 
we remove a large number of broad-line AGNs below our detec- 
tion threshold. However, even when we include these objects, 
our luminosity function flattens at considerably higher luminos- 
ity than the Hao et al. function. Relative to Hao et al. we have 
removed objects both with our initial Ha detection algorithm 
and with our line-fitting procedure, at the <~ 20% level. Be- 
cause the luminosities involved are very low, the corresponding 
volume corrections are very large, leading to highly discrepant 
space densities in the two studies. In some respects it is not 
surprising to find such a large range of possible space densities 
at these very low Ha luminosities. As we show in Appendix 
A, the measurements of Lu a are very inaccurate in this regime. 
Perhaps the most important message of this exercise is that the 
SDSS is not an adequate tool to measure the line luminosity 
function below <~ 10 40 5 ergs s" 1 . 

In addition to the differences in object selection, we adopt 
a very different set of assumptions in building our luminos- 
ity functions than Hao et al. They have used the maximum 
likelihood formalism, but in calculating their selection function 
[pi(L)] have assumed that neither the line shape nor the galaxy 
luminosity are linked to the Ha luminosity. In our V /V ma x for- 
malism, we need not make any assumption about galaxy lu- 
minosity. However, when we calculate our selection function 
for the maximum likelihood function, we do not allow a given 
object to exceed its Eddington luminosity, thus implicitly as- 
suming a correlation between galaxy and AGN luminosity. For 
the majority of objects, we are thus adopting a higher maxi- 
mum luminosity than the fiducial luminosity of ~ 10 42 ' 6 ergs 
s" 1 assumed by Hao et al.. Furthermore, we assume that as the 
AGN luminosity decreases for a given object, the line width in- 
creases, thus increasing the noise per pixel across the broad line, 
and decreasing our detection efficiency at low luminosity. The 
net effect at low luminosity is apparently that we have larger 
volume corrections than Hao et al., leading to lower inferred 
space densities. 

As an aside, we note that our luminosity function is consis- 
tent with that of Croom et al. (2004), with the break luminosity 
~ 1 mag lower. On the other hand, while the Hao et al. function 
appears to be consistent with that of Ulvestad & Ho (2001; see 
Fig. 1 1 of Hao et al.), the galaxy luminosity is included in the 
latter luminosity function. In most cases the AGN luminosity 
itself is significantly fainter, and it is probably fair to shift their 
function fainter by ~ 4 magnitudes (see Ho 2004 for details). 

4.3. X-ray Luminosity Function 

As an additional sanity check, we compare our luminosity 
function with the zero-redshift luminosity functions presented 
by Hasinger et al. (2005). This work focuses predominantly 
on soft X-ray-selected samples with broad Balmer lines in op- 
tical follow-up (see also Schmidt et al. 1998). To perform a 
comparison, we (statistically) convert our observed Ha lumi- 
nosities to soft X-ray luminosities. We employ the Greene & 
Ho (2005b) conversion between Ha luminosity and L5100, and 
then, assuming a spectral slope of a = 0.44 (f v oc v~ a \ see, e.g., 
Vanden Berk et al. 2001; Greene & Ho 2005b), the relation 
between a ox and L2500 from Strateva et al. (2005). Although 
the Strateva et al. relation (see also Steffen et al. 2006) is only 
measured for UV luminosities of ~ 10 43 — 10 48 ergs s" 1 , Greene 
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& Ho (2007) find that, in general, the low-mass AGNs from 
Greene & Ho (2004) obey the extrapolation of the Strateva re- 
lation to lower luminosity. Because there is large scatter in the 
conversion between X-ray and optical luminosity, we generate 
1000 X-ray luminosity functions from our optical function, in 
each case perturbing the intercept and slope of the relation by a 
log-normal deviate determined by the measured scatter in each 
variable. We then select the median value at each luminosity, 
and the upper and lower bounds to bracket 68%. The resulting 
inferred soft X-ray luminosity function is shown in Figure 6, 
with the Hasinger et al. double power-law fit shown as a solid 
line. Given the uncertainties involved, we find quite reasonable 
agreement. 

Table 3. BH Mass Function 



log (M BH ) 


log (i) 


(1) 


(2) 



5.00 
5.25 
5.50 
5.75 
6.00 
6.25 
6.50 
6.75 
7.00 
7.25 
7.50 
7.75 
8.00 
8.25 
8.50 
8.75 



a 1 n + 0-16 

-6.10T 25 

-5 81+° 17 

-6 05+° 11 
D - uo _0.14 

-5 51+ - 07 
°- ol _0.08 
- 91 +0.08 

-5.06+ 010 



-0.13 
3+0.05 
-0.05 
-5 10+ 004 

-5 10+ 005 
°- iu -0.05 
, +0.03 
-0.03 

-5 50+ ' 06 
°- ou -0.07 

-5 79+ 06 
J ' ' -0.07 
,+0.05 
-0.06 
R OO+0.05 

,+0.09 
-0.11 
3+0.11 
-0.14 



-5.13" 1 



-5.31 1 



-6.22" 1 



-7.32] 
-7.98" 1 



Note. — Col. (1): Virial 
Mbh (Mq) as inferred from 
L H « and FWHM H „. Col. 
(2): Space density (# Mpc -3 
log Mbh -1 ) as a function of 



4.4. BH Mass Function 

Using the V mdx weights derived above, we show in Figure 7 
volume-weighted BH mass function (the number of BHs per 
unit volume, per unit logarithmic BH mass, in bins of 0.25 
dex in BH mass; see also Table 3). As with the luminosity 
function, we fit the mass function with a Schechter function, a 
double power law, and a log-normal parameterization. In this 
case we include all mass bins in the fitting. Once again, the 
Schechter function provides a poor fit, with \ 2 r = 18- Both 
the double power-law and log-normal functions provide rea- 
sonable fits, with xf — 2-5 and xf = 1-2, respectively. The 
best-fit double power law is shown as a dashed line in Figure 
7, and is parameterized by </>(M») = (4.7 ± 0.5) x 10" 6 Mpc" 3 , 
M* = 10 732±005 Mq, a = 0.78 ± 0.06, and /3 = 3.00 ± 0.07. 
The log-normal fit, shown as a solid line in Figure 7, has a 
normalization of (6.4 ±0.3) x 10~ 6 Mpc" 3 , a central mass of 
10 6.6±o.03 Mq ^ and a wMth f o.59 ±0.01 dex. 

Taken at face value, our BH mass function displays a clear 
break at ~ 10 6 6 M , and declines toward both higher and lower 
mass. This is comparable to the characteristic mass of local 



narrow-line AGNs (Heckman et al. 2004) but is significantly 
lower than that of inactive BHs (~ 10 8 M Q ; e.g., Marconi et al. 
2004). Unfortunately, as demonstrated by Figure 3, our data are 
dominated by the strong selection bias inherent in a magnitude- 
limited sample. It is clear that, in terms of luminosity, we are 
significantly incomplete below ~ 10 40 75 ergs s" 1 . However, in- 
completeness in luminosity is redshift-dependent and may not 
be translated directly into incompleteness as a function of BH 
mass. At a given L Hq , the S/N across the line (and thus its de- 
tectability) decreases with increasing FWHMh q , and thus in- 
creasing BH mass. Additionally, as the BH mass increases, the 
corresponding bulge luminosity is presumably increasing, fur- 
ther decreasing the contrast and the S/N in the broad line. 
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Fig. 7. 
log M B h 



log (M BH / M ) 

- Volume-weighted BH mass function in bins of 0.25 dex (# Mpc -3 
1 ). The weights used are identical to those for the luminosity func- 
tion, and as above we show in the inset the mass functions for objects targeted 
as galaxies (left) and QSOs (right), respectively. Although we are subject to 
significant incompleteness, we will argue below that there is truly a turnover 
in active galaxy masses at both lower and higher BH masses. We have fit the 
mass function with both a double power-law (dashed line) and a log-normal 
function (solid line). 

If we could uniquely ascribe a host galaxy luminosity and 
light profile to each Mbh, then we could easily model our in- 
completeness as a function of M BH and L\, i/L E ^. At high BH 
mass, this is in fact possible, as there is a relation linking Mbh 
and spheroid luminosity, and the fundamental plane tells us the 
typical sizes (and thus fiber luminosities) of elliptical galax- 
ies. However, for spiral or dwarf spheroidal host galaxies, there 
ceases to be a unique mapping between M B h and galaxy lumi- 
nosity or structure. The relation between bulge-to-total ratio 
and galaxy luminosity is poorly quantified and contains sig- 
nificant scatter in any case. Furthermore, at lower masses, as 
the AGNs become intrinsically fainter, only systems with rel- 
atively luminous host galaxies will fall above the magnitude 
limit of the SDSS. For these reasons, at low M B h the calculated 
incompleteness is a strong function of the assumed (but uncon- 
strained) host galaxy morphology. 

As a matter of practicality, then, we turn the problem around. 
Rather than attempting to quantify our absolute incomplete- 
ness as a function of Mbh, we simply quantify the range of 
host galaxy luminosities for which we might hope to detect a 
BH of a given mass, L Ha , and z. Over narrow ranges in all of 
these parameters, neither the line width nor the galaxy contin- 
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Incompleteness as a function of galaxy luminosity for simulated galaxies. These simulations are run for BH masses of Mbh= 10 



■ 10 6 4 M© at three 



redshifts (z = 0.04, 0.07, 0.10), and with the limited range of L\\ a = 10 40,5 — 10 41 ergs s _1 . Although the input fiber luminosities are — 14 < Mg < -22, there are 
non-zero detection fractions only over the range of fiber luminosities shown here. The simulations demonstrate that we have roughly constant incompleteness over 
this range in Mbh, ^-Hct, and I- 



uum strength changes dramatically and the completeness ought 
to be constant (provided the BHs are drawn from the same host 
galaxy population). In such bins, with uniform completeness, 
we are able to measure true changes in space density. Simula- 
tions allow us to isolate ranges of M BH , L Hq , and z with constant 
sensitivity to galaxy fiber luminosity. Note, however, that in 
any given interval, we necessarily exclude different members 
of the population as a function of Mbh; at the lowest masses 
we preferentially exclude those systems in faint hosts, while 
at the highest Mbh we exclude the higher Lboi/^Edd systems. 
Therefore, we implicitly assume that the distribution of Mbh 
is uniform independent of both disk luminosity and Eddington 
luminosity. With this approach, however, we need not concern 
ourselves directly with host galaxy structural per se, but solely 
the luminosity. The total host galaxy luminosity must be high 
enough that the source is spectroscopically targeted, while the 
fiber luminosity must be low enough to allow detection of the 
broad line. 

We investigate three mass regimes (Mbh=10 5 5 - 10 6,4 , 10 6 - 
10 7 , and 10 6 5 - 10 7 5 M@), choosing optimal z and Lhq ranges 
for each. Our procedure is described in most detail for the low- 
est (and most challenging) mass bin, and then results are pre- 
sented for all three. In the first bin, we are fundamentally lim- 
ited by the total number of objects. Therefore, we are forced 
to use the lowest possible redshift bins: z = 0.05 - 0.07 and 
z = 0.07-0.1. These are bins with > 10 objects per bin for 
the most part; at still lower z a prohibitively large range in 
distance is needed to populate each bin. In terms of Ln a , the 
highest luminosity is set by the Eddington luminosity of the 
lowest mass bin, in this case 10 41 ergs s" 1 for a BH with mass 
10 5 5 M Q , while the lowest luminosity is set by the paucity of 
lower-luminosity objects (10 40,5 ergs s -1 ). 

Simulations allow us to verify that the selection probabil- 
ity is indeed independent of host galaxy fiber luminosity. We 
make artificial spectra in the appropriate Mbh and Lhq range, 
with fiber galaxy luminosities spanning -14 < Mb < -22. The 
galaxy continuum is modeled as a single stellar absorption-line 
system, constructed from the eigenspectra of Yip et al. (2004), 
and the S/N is varied to correspond to typical SDSS spectra 

3 At a given luminosity, a wide range of galaxy morphologies are permitted. Thus it is still possible to find significant differences in <r* (e.g., Greene & Ho 2006b) 
and potentially host galaxy structure (J. E. Greene, in preparation) as a function of Mbh for the SDSS-selected samples of low-mass systems. 



over the redshift range of interest. Five realizations are made 
for each galaxy luminosity and S/N, and each spectrum is run 
through our full detection algorithm. For those with detectable 
broad Ha, we then investigate whether the galaxy luminosity is 
sufficient for spectroscopic targeting in the first place. Recall 
that this limit depends on total (rather than fiber) luminosity, 
but there is not a one-to-one conversion from fiber to total lu- 
minosity; it depends on galaxy morphology and redshift rather 
strongly. Therefore, we place an upper limit on the total galaxy 
luminosity by insisting that the fiber luminosity account for no 
less than 20% of the total galaxy luminosity (as motivated by 
the observed range shown in Fig. 9 of Tremonti et al. 2004). 
Over the entire range of galaxy luminosities we explore, a non- 
zero detection fraction results only for fiber luminosities in the 
range -16 < Mb < -18, but the detection fractions at a given 
host luminosity are very constant across the mass range of in- 
terest, as shown in Figure 8. 

In Figure 9a we show the resulting mass functions for the 
two different redshift bins. Visually, it appears that the space 
density is truly falling at low mass. To quantify the significance 
of this result, we perform a simultaneous fit to the mass func- 
tions in the two redshift bins. We have chosen to use a single 
power-law in this case, 

0*(M„) / M 
^ M* \M* 

The two redshift bins are constrained to have the same value of 
a but are allowed different values of </>*(M»). An apparent flat 
slope in our logarithmic mass function, $(Mbh), would corre- 
spond to a =-1. We find a reduced x?=101 and a best-fit slope 
of a = -0.6 ± 0.2, corresponding to a 2 a detection of a falling 
space density for M BH < 10 6 5 M Q . If we relax the luminos- 
ity criterion to include all objects with Ln a = 10 40 - 10 41 ergs 
s , the best-fit slope is a = -0.7 ±0.1. We have performed this 
fit including the lowest-mass point (M B h= 10 5 ' 2 - 10 5 ' 5 M Q ), 
but this point is suspect both because there only 5 objects in 
the bin, and also because our maximum Ln a exceeds the Ed- 
dington limit for BHs in this bin. However, when this point is 
removed, the best-fit slope does not change. The major out- 
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standing caveats, unavoidable products of our technique, are 
that these BHs must both be highly active and live in fairly lu- 
minous host galaxies. 

We repeat the above analysis for two other mass intervals. 
For the first, Mbh = 10 6 — 10 7 Mq, we employ a luminosity in- 
terval of Lua = 10 40,5 - 10 417 ergs s" 1 and, because there are 
many more objects in this mass range, three redshift intervals, 
z = 0.05 - 0.08, 0.08 - 0. 1 1 , and 0. 1 1 - 0. 14. The fitting results 
are shown in Figure 9b, and, with a best-fit a = -0.72 ±0.06, are 
consistent with the slope found above. In this case, we have a 
5 a detection of decreasing space density below M Bli < 10 7 M©. 
However, it should be noted that the best fit has x 2 r = 2.7, sug- 
gesting that this is not a terribly good model. While we might 
achieve an improved fit with a broken power-law model, this 
hardly seems justified for such a limited range in BH mass. 
Finally, we turn to the break at high mass, using the mass 
bin M B h = 10 6 5 - 10 7 ' 5 Mq, and the luminosity bin L na = 
10 41 - 10 42 ergs s" 1 . Again we are able to utilize three redshift 
bins. Because of the bright magnitude limit of the SDSS, we 
are forced to go to higher redshift for completeness in the high- 
est mass bins here, z = 0.10-0.13, 0.13-0.16, and 0.16-0.19. 
In this mass bin, the mass function begins to turn over: we find 
a best-fit a = -1.13 ±0.04, corresponding to a 3 a detection 
of a falling space density above M Bii > 10 6 5 Mq. Again, the 

single power-law model is not a fantastic fit (x?=3.4), so it is 
difficult to say when the turnover actually occurs. If we repeat 
the fit without the highest mass bin, the fit apparently improves 
(Xr = 1-6) and the slope steepens, a = -1.40 ±0.05. Even if 
we neglect the top two or three mass bins, we find the same re- 
sult. Therefore, we cannot exactly localize the break; it occurs 
somewhere in the range 10 6 5 — 10 7 5 Mq. 



4.5. Comparison with Heckman et al. (2004) 

For the first time, we are able to compare narrow- and broad- 
line AGN mass functions rather than just luminosity functions 
(Fig. 10). The narrow-line objects come from the SDSS 
database of Kauffmann et al. (2003), for which BH masses 
are inferred from the stellar velocity dispersion (Heckman et 
al. 2004; we actually use the most recent sample from DR4). 
Reliable BH masses extend no lower than 10 6 ' 3 Mq due to the 
resolution limit of the SDSS (a ~ 71 km s" 1 ; e.g., Heckman et 
al. 2004). Volume weights depend only on galaxy luminosity, 
since they were selected from the main galaxy sample. 

It is clear that the space density of narrow-line AGNs in the 
Heckman et al. sample is higher than that of the broad-line sam- 
ple by an order of magnitude. However, it is crucial to remem- 
ber that the selection effects for the two samples are radically 
different. In particular, narrow emission lines may be detected 
to significantly lower luminosities than broad lines. Therefore, 
while the typical Eddington ratio in our sample is Lboi/^Edd 
« 0.1, the typical L bol /L Edd for a 10 7 Mq BH in the Heck- 
man et al. sample is closer to Lboi/£Edd~ 0.04 (see their Fig. 
3). It would be much more informative to compare the broad- 
and narrow-line objects in matching mass and luminosity bins. 
However, to do so properly requires a consistent measure of 
AGN luminosity across the two samples. Typically, the [O III] 
line is used (e.g., Zakamska et al. 2003), but for broad-line ob- 
jects there may be significant contamination to this line from 
broad Fe II emission (e.g., Greene & Ho 2005b). Since we have 
not properly subtracted the Fe II, a more detailed comparison 
must await future work. 



5. DISCUSSION 
5.1. Active Fraction or Duty Cycle 

We have constructed a BH mass function for the local Uni- 
verse, using broad-line AGNs and virial BH mass estimates. 
By looking at very narrow ranges in redshift and AGN lumi- 
nosity, we have detected a break at M BH « 10 6 5 - 10 7 Mq. We 
now investigate what fraction of BHs are active as a function 
of BH (spheroid) mass, by computing the inferred shape of the 
inactive BH mass function, using the formalism of Marconi et 
al. (2004; see also Yu & Tremaine 2002; Shankar et al. 2004; 
McLure & Dunlop 2004). Statistically speaking, the fraction of 
currently active systems may be interpreted as a duty cycle for 
BHs in that mass range. 

Marconi et al. compute consistent inactive mass function 
from a variety of optical and near-infrared surveys and so for 
simplicity we use the Kochanek et al. (2001) /f-band galaxy 
luminosity function. The galaxy luminosity and BH mass func- 
tions are linked through the /T-band M B h - ^buige relation of 
Marconi & Hunt (2003), assuming an intrinsic scatter in that 
relation of 0.3 dex and using Monte Carlo realizations to incor- 
porate measurement uncertainties. The resulting uncertainties 
(see their Fig. 2b) are relatively constant for M B n< 10 8 Mq, 
and have an amplitude of a « 0.2 dex, although there are a va- 
riety of systematic uncertainties that preferentially impact the 
low-mass regime (for comparison, the Shankar et al. mass func- 
tion is ~ 0.5 dex higher than Marconi et al. at M B h= 10 7 M q ). 
At high mass, the mass function is well-constrained; the lumi- 
nosity of elliptical galaxies may be translated directly into BH 
mass. On the other hand, to properly translate between galaxy 
luminosity and BH mass at low luminosities requires knowl- 
edge of the (poorly constrained) luminosity function of bulges. 
For the moment, we simply adopt the fiducial values from Mar- 
coni et al., who assume morphological fractions (mb u i ge - m tot . d i) 
fromFukugitaet al. (1998) of 1 1% elliptical, 21% SO, 43% Sab, 
and 19% Scd and bulge-to-disk values from Aller & Richstone 
(2002; derived from the data of Simien & de Vaucouleurs 1986) 
of 0.64 ± 0.30 for SO, 1 .46 ± 0.56 for Sab, and 2.86 ± 0.59 for 
Scd galaxies. We show the resultant distribution of inactive BH 
mass as a dashed line in Figure 11a. For reference, the con- 
tribution from early- and late-type spirals within the adopted 
formalism are shown in long and short dot-dash lines, respec- 
tively. 

The active fraction clearly depends strongly on BH mass. 
The differences in shape between active and inactive systems 
is highlighted in Figure 11a, where the active mass function is 
boosted by a factor of ~ 200 to overlay the inactive mass func- 
tion derived as described above. (Because the low-mass slope 
is particularly uncertain, we also show the best-fit slope to the 
lowest-mass and redshift bin from §4.4 [dotted]). In the bot- 
tom panel we show the ratio of each subpopulation with the 
inactive mass function. BHs with masses ~ 10 7 Mq, living in 
~ 10 10 M Q spheroids are active ~ 0.4% of the time. At higher 
BH mass the number of active BHs drops quickly; local mas- 
sive BHs are not radiating at substantial fractions of their Ed- 
dington limit. The same behavior was seen by Heckman et al. 
(2004) for narrow-line AGNs. It is reassuring that we find con- 
sistency with the results of Heckman et al., despite our differing 
(indirect) means of measuring Mbh- The Heckman et al. result, 
and now ours, may be seen as confirmation of the suggestion, 
based on AGN luminosity functions, that lower-mass objects 
build their mass at low redshift, while the most massive BHs 
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FIG. 9. — Fits to the BH mass function in narrow ranges of Mbh, ^Ha, and z, as indicated. These windows were selected to have uniform sensitivity to galaxy 
fiber luminosity, so that the measured slopes represent physical changes in space density. We perform simultaneous power-law fits (allowing different amplitudes but 
fixing the slopes between bins) to the redshift bins shown in each panel. We clearly detect a break in the mass function in the mass range 10 6 5 - 10 7 Mq, although 
the exact location is not clear from these data. Note also that we are only sensitive to active galaxies with a narrow range of fiber luminosities, which may bias these 
results (see text). 



completed their growth at high redshift ("cosmic downsizing"; 
Cowie et al. 1996; Ueda et al. 2003; Steffen et al. 2003). 

There is an important caveat that deserves discussion. While 
Marconi et al. (as well as Shankar et al. and McLure & Dunlop) 
find good agreement between inactive mass functions derived 
from galaxy luminosity and er* functions, Tundo et al. (2007) 
report that the mass function based on er» distributions is sys- 
tematically lower than that based on galaxy luminosity func- 
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tions. 

FIG. 10. — Comparison between the BH mass function of broad-line (filled 
symbols; this paper, truncated at 10 6 Mq for consistency with the narrow- 
line AGN sample) and narrow-line (dashed line; Heckman et al. 2004) AGNs. 
Mass function for broad-line objects as described above. Volume corrections 
for the narrow-line objects are based solely on the galaxy color and luminosity. 
The space density of narrow-line objects in their sample is significantly higher 
than what is seen for broad-line objects, due to the relative ease of detecting 
high-EW narrow lines. 



According to Bernardi et al. (2007), the M B h - c* calibration 
sample (e.g., Tremaine et al. 2002) has a biased Faber-Jackson 
(1976) relation compared to the SDSS, in the sense that the 



galaxies have lower luminosities at a given 0% (or higher a* at 
a given L). One consequently infers a lower BH mass function 
from the M Bll -a* relation than from the M B H-£buige relation. 
Bernardi et al. (2007) present a model in which the Mbh-c* re- 
lation provides the more reliable estimator of BH mass density. 
If (and this is by no means certain) they have properly mod- 
eled the bias, then the true inactive BH mass density should be 
~ 0.5 dex lower at M B h= 10 9 Mq, and the active fraction de- 
clines less steeply with mass than we report here. At the present 
time, given the systematic differences between authors and es- 
timators, we conservatively assign a factor of 3 uncertainty to 
our active fraction across all mass bins. 

It is also worth noting that the active fraction depends sen- 
sitively on the depth of the survey. The Palomar spectroscopic 
survey of galaxies, which is significantly more sensitive to low- 
luminosity Ha lines, finds an active fraction of ~ 60 % in 
nearby, bulge-dominated galaxies (Ho et al. 1997b). Despite 
the limited dynamic range in our sample, the number density of 
sources at a given BH mass increases at lower Eddington ratio, 
modulo incompleteness (Fig. lib). 

Aside from the known decrease of active fraction at high BH 
mass, we additionally see a decline in active fraction for BHs 
< jo 6 - 5 M©, but unlike the case of massive BHs, we cannot 
clearly interpret the decline. Taken at face value, the observa- 
tion suggests that as we move to masses below the peak, the 
duty cycle declines. Unfortunately, both the active and inactive 
mass functions are highly uncertain for masses < 10 6,5 Mq. Ul- 
timately, we are limited by our knowledge of (a) the shape of 
the luminosity function, (b) the conversion between galaxy lu- 
minosity and BH mass for systems that contain BHs, and (c) the 
occupation fraction of BHs in dwarf galaxies. We discuss the 
potential magnitude and impact of each of these effects below. 

First of all, the luminosity function of Kochanek et al. (2001) 
is measured only to a luminosity of M K = -20 mag, or a BH 
mass of 10 6 4 Mq for an elliptical galaxy, or 10 5 ' 1 Mq for an 
Scd galaxy. Thus, the entire low-mass regime involves some 
extrapolation. Deeper luminosity functions do exist (e.g., Blan- 
ton et al. 2005), but we are still left with the highly uncertain 
task of converting between an observed luminosity function and 
a BH mass function. The morphological fractions adopted by 
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Marconi et al. (2004; from Fukugita et al. 1998) are relatively 
well-determined for E/SO galaxies. They are adopted from the 
morphology-density study of Postman & Geller (1984), based 
on the CfA redshift survey (Huchra et al. 1983; morphologies 
from de Vaucouleurs et al. 1976 and Nilson 1973), which sep- 
arated galaxies into E, SO, and spiral using visual classifica- 
tion. On the other hand, the fractions of Sab and Scd galaxies 
are derived in a less-direct fashion. Tinsley (1980) took the 
observed color and apparent magnitude distributions from the 
photographic plate study of Kirshneret al. (1978) and solved si- 
multaneously for the morphological fractions and redshifts that 
reproduced the observations, using broad-band color distribu- 
tions from Pence (1976). Because of the degeneracies inherent 
in photometric redshifts based on two broad-band colorfl and 
the limited sample size (~ 800 galaxies), this leads to signifi- 
cant uncertainties in the bulge fractions for spiral galaxies. 

Ideally, we would like a direct measurement of bulge (rather 
than total galaxy) luminosity. Dzanovic et al. (2007; see also 
Tasca & White 2007) present spheroid luminosity and mass 
functions based on two-dimensional image decomposition of 
~ 9000 SDSS galaxies. The resulting BH mass function is con- 
sistent with a decrease in BH mass density of ~ 0.5 dex be- 
tween M B h = 10 7 M and M B h = 10 5 M Q . However, the con- 
version from spheroid mass to BH mass relies on extrapolation 
of BH-bulge relations that have been calibrated only for higher 
masses. The host galaxies of low-mass BHs (Mbh < 10 6 Mq) 
may not even contain bulges. Indeed, neither of the two best- 
studied low-mass BHs (M BH « 10 5 Mq) — NGC 4395 (Filip- 
penko & Ho 2003) or POX 52 (Barth et al. 2004)— contains 
a classical bulge. The former is a late-type, bulgeless spiral, 
which lies on the low-mass extrapolation of the MBH-^tmige re- 
lation only if one considers the luminosity of its nuclear star 
cluster. POX 52 is a dwarf spheroidal galaxjO, and from the to- 
tal luminosity of the host one would predict a BH mass 50 times 
higher than its virial mass or cr* suggests. A similar trend is 
seen for the SDSS sample of Greene & Ho (2004; J. E. Greene 
et al., in preparation). It is possible that these more distant sys- 
tems also obey an MBH-^buige relation with a nuclear star clus- 
ter, and POX 52 may show evidence of nucleation from HST 
observations (C. J. Thornton et al., in preparation), but in gen- 
eral it is not possible to predict nuclear cluster luminosity from 
galaxy luminosity (although the two may be correlated; Fer- 
rarese et al. 2006; Rossa et al. 2006; Wehner & Harris 2006). 
In short, a conversion from galaxy luminosity to BH mass is 
presently unconstrained for BHs with masses Mbh< 10 6 Mq. 

Finally, while we believe that most (if not all) bulges host 
BHs, the occupation fraction of BHs in late-type galaxies in 
unknown. Dynamical measurements in the Local Group for 
the late-type galaxy M33 (Gebhardt et al. 2001) and the dwarf 
spheroidal galaxy NGC 205 (Valluri et al. 2005) show that nei- 
ther hosts a BH with a mass in accordance with expectations 
from the Mbh - cr* relation. On the other hand, the massive 
globular cluster Gl in Andromeda hosts a BH in accord with 
the M BH - o* relation (Gebhardt et al. 2002, 2005), and NGC 
4395, POX 52, and the Greene & Ho (2004) sample, all late- 
type galaxies hosting BHs with masses < 10 6 Mq, are consis- 
tent with the low-mass extrapolation of the Mbh _ cr* relation 
(Filippenko & Ho 2003; Barth et al. 2004, 2005). 



In light of these many uncertainties, it is difficult to inter- 
pret the apparent decrease in active fraction at low mass. The 
most straightforward interpretation may be that the shape of the 
active mass function actually traces the shape of the bulge lumi- 
nosity function at low luminosity. Otherwise, is it is difficult to 
understand why late-type spiral galaxies, with a ready gas sup- 
ply, nevertheless are preferentially inactive. The distribution of 
activity in the Palomar spectroscopic survey of nearby galax- 
ies support this hypothesis; Ho et al. (1997b) find that Seyfert 
galaxies (i.e., objects with relatively high accretion rates) are 
predominantly found in early-type spiral galaxies (Sa-Sb). In 
the future, with the availability of more complete measurements 
of spheroid luminosity functions, we may use our observations 
of active fractions to measure the mass at which the BH occu- 
pation fraction departs from unity. 

5.2. Connection to Higher Redshift 

Broad-line AGN activity at the present day is dominated by 
~ 10 7 Mq BHs radiating at — 10% of their Eddington lim- 
its (Fig. Mb). During the quasar epoch, 1 < z < 3, BH 
mass growth is dominated by optically bright, near-Eddington 
growth of massive BHs (e.g., Vestergaard 2002) and this growth 
basically accounts for most of the total BH mass density ob- 
served today (e.g., Yu & Tremaine 2002). Even as one moves 
to lower luminosities, BH growth appears to be dominated by 
near-Eddington accretion, at least for massive BHs at high red- 
shift (Kollmeier et al. 2006). In general, this picture sits com- 
fortably with the observation that massive elliptical galaxies 
(hosts of the most massive BHs) formed their stars rapidly 
at high redshifts (e.g., Bower et al. 1992; Trager et al. 2000; 
Thomas et al. 2005). It is the redshift range between the present 
and z«l that remains poorly constrained. From the point of 
view of galaxy evolution, there appears to be a significant (fac- 
tor of 2) increase in the mass density of spheroids over this 
interval (e.g., Brown et al. 2007; Faber et al. 2007). It would 
be interesting to know whether any commensurate BH growth 
occurred. 

There are tantalizing hints. Optical QSO luminosity func- 
tions show that massive BHs are not radiating at high Eddington 
ratio below z ~ 1 (e.g., Richards et al. 2006), but do not con- 
strain the evolution of low-luminosity sources. X-ray luminos- 
ity functions, built from very deep Chandra and XMM-Newton 
pointings, are able to probe significantly fainter luminosities 
(see Brandt & Hasinger 2005 for a recent review). It is clear 
that the space density of low-luminosity (Lx ~ 10 42 - 10 43 ergs 
s" 1 ), X-ray-selected sources peaks at considerably lower red- 
shift (z k 0.7) than higher-luminosity sources. Unfortunately, 
multiple scenarios for BH growth result in similar luminosity 
evolution. At one extreme, one may imagine all near-Eddington 
growth at all masses occurs at high redshift and then slowly 
shuts off, while at the other extreme, one may imagine that the 
low-luminosity radiation is dominated by the growth of low- 
mass BHs, which acquire a significant fraction of their mass at 
low redshift. In reality, some combination of these two scenar- 
ios presumably occurs. The latter hypothesis, dubbed "cosmic 
downsizing," is often preferred in the literature, partially be- 
cause galaxies display similar behavior (e.g., Cowie et al. 1996; 
Barger et al. 2005). 

4 For instance, Fukugita et al. (1995) find that the X'-corrections for Sa-Sc galaxies at z < 0.3 differ by Sz < 0.1, in B or V. Degeneracies between metallicity and 
age only add to the uncertainties. 

5 Such objects are commonly referred to as dwarf elliptical galaxies in the literature. However, because their structure is quite different from that of classical 
elliptical galaxies, we prefer to refer to them as dwarf spheroidal systems. 
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FIG. 11. — (a) Comparison of the inactive BH mass function with the broad-line BH mass function. The inactive mass function is derived from the galaxy 
luminosity function of Kochanek et al. (2001), using the formalism of Marconi et al. (2004), including an intrinsic dispersion of 0.3 dex in the relation between 
galaxy luminosity and BH mass. We show the log-normal fit to the active BH mass function offset by 2.35 dex (solid) and the best-fit power law to the lowest-mass, 
lowest-z bin (Fig. 9a) to demonstrate the uncertainty in slope and break position at the low-mass end. The ratio between the total inactive mass function and each 
population is shown in the bottom panel. The vertical line indicates the point at which we begin to use an extrapolation of the luminosity function, (b) Comparison 
of the inactive BH mass function (dashed) with the active mass function in different Eddington ratio bins (solid symbols). Bins with > 3 points are shown. It is clear 
that all high-mass BHs are radiating at very low Eddington ratios. 



The first concrete evidence in support of downsizing came 
from the Heckman et al. (2004) study, which showed unam- 
biguously that local accreting BHs are typically an order of 
magnitude less massive than the typical local inactive BH. The 
present study finds the same result. However, we still must ask 
whether X-ray-selected AGNs at z = 0.7 are really the same 
population as local optically selected sources. The story may 
not be so straightforward. Netzer & Trakhtenbrot (2007) find 
that the typical BH mass and Eddington ratio of optically se- 
lected AGNs both increase significantly from the present time 
to z = 0.7, so that what we measure locally may no longer apply 
at that epoch. The other intriguing, but not definitive, evidence 
comes from recent optical studies of the host galaxies of X-ray- 
selected AGNs. It seems that the majority of sources are found 
in massive, relatively red galaxies (e.g., Colbert et al. 2005; 
Nandra et al. 2007; see also Barger et al. 2005). The clustering 
properties of X-ray-selected AGNs also suggest that the hosts 
are typically more massive than optically-selected AGNs at the 
same redshift (e.g., Coil et al. 2007, Miyaji et al. 2007, Coil, 
A. L. private communication). Assuming that the M Bll -<j* re- 
lation is largely in place at z ~ 0.7 (Peng et al. 2006; Salviander 
et al. 2007; but see also Treu et al. 2004; Woo et al. 2006), the 
observed galaxy luminosities correspond to BHs with masses 
> 10 8 Mq, casting doubt on the cosmic downsizing scenario 
for this population. 

Taken at face value, the fact that X-ray-selected samples 
are found in very luminous host galaxies suggests that low- 
luminosity, hard X-ray-selected sources consist predominantly 
of massive BHs in a low-Eddington state. These sources are 
known to be an intrinsically hard population with no big blue 
bump (Ho 1999), and if they have broad lines, they will be ex- 
tremely low contrast, and thus difficult to detect. The conjec- 
ture that low-Eddington sources dominate the hard X-ray popu- 
lation thus accounts for both the general lack of apparent broad 



emission and the generic hard spectral shape of the X-ray back- 
ground (see additional arguments in Shen et al. 2007). On the 
other hand, low-luminosity, soft X-ray-selected sources, most 
likely closely linked to optically selected broad-line AGNs, also 
have a peak space density at low redshift (e.g., Hasinger et 
al. 2005). At the present time, it is unclear whether the mas- 
sive host galaxies measured in the studies above are represen- 
tative of the population as a whole. An unbiased comparison 
of the Eddington ratio distributions of the soft and hard X-ray- 
selected samples is needed. 

A rather different approach to constraining the evolution of 
BH mass density is taken by Merloni (2004). He constructs a 
BH mass function using constraints from the joint X-ray and 
radio luminosity function, in combination with an empirical re- 
lation between BH mass and X-ray and radio luminosity (the 
"fundamental plane"; e.g., Merloni et al. 2003). He finds fast 
evolution; by z = 0.6 the characteristic accreting BH mass has 
increased from 10 7 M Q to 10 8 Mq. However, his analysis also 
tends to favor low-Eddington sources, and indeed his charac- 
teristic L x /L E dd never exceeds 1%. Ideally, work like that pre- 
sented here may be extended to higher redshift to directly ad- 
dress the distribution of BH mass at intermediate redshift. 

6. SUMMARY 

We present the first measurement of the local BH mass func- 
tion for broad-line active galaxies. Our sample of <~ 9000 galax- 
ies is drawn from the Fourth Data Release of the Sloan Digital 
Sky Survey, based on the presence of broad Ha emission. Our 
Ha luminosity function is found to be consistent with previous 
measurements, and with the local soft X-ray luminosity func- 
tion. Using standard scaling relations between AGN luminos- 
ity and line width, we derive BH masses for the entire sam- 
ple. Much like the Heckman et al. (2004) sample of narrow- 
line AGNs, the typical BH in our sample has a mass of M B h 
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« 10 7 Mq and an Eddington ratio of Lboi/^Edd ~ 10%. Al- 
though we are highly incomplete for low-Eddington ratio sys- 
tems, by looking at very narrow bins in Mbh, £hq, and z we 
find clear evidence for a true turnover in space density below 
Mbh ~ 10 6 5 - 10 7 Mq. Compared to the inferred shape of the 
mass function of inactive BHs, the mass function of active BHs 
falls significantly both above and below this characteristic break 
mass. The dearth of active massive BHs is a familiar result — 
massive BHs are mostly quiescent in the local Universe. The 
decreasing space density at low BH mass presumably reflects 
the fact that bulge fraction and BH occupation fraction both de- 
crease in dwarf galaxies. 

One of the major motivations of this study was to con- 
strain the shape of the BH mass function for BHs with masses 
< 10 6 Mq, since broad-line AGN surveys currently have the 
unique capability to explore this mass regime. We have found 
evidence for a decreasing space density of active BHs with 
masses < 10 6 5 Mq, at least for objects with host galaxies more 
luminous than Mb ~ -16 mag. However, substantial work re- 
mains to determine the true demographics of BHs in low-mass 
galaxies. Measurements that isolate the luminosity and mass 
functions of spheroids (e.g., Dzanovic et al. 2007) on the one 
hand, combined with empirical conversions between spheroid 
luminosity and BH mass in the low-mass regime (J. E. Greene 
et al., in preparation) on the other, will provide predictions for 
the shape of the inactive BH mass function for Mbh < 10 Mq. 
At the same time, alternate search techniques, which are less 
sensitive to host galaxy luminosity, are required to discover 



low-mass, low Lboi/^Edd BHs. Apart from deeper nuclear op- 
tical spectroscopic surveys, deep X-ray and radio surveys may 
provide sensitivity to lower-Eddington ratio systems and elimi- 
nate the host luminosity bias suffered by the present study. Fi- 
nally, future time-domain surveys may place limits on the fre- 
quency of flares from the tidal disruption of stars captured by 
BHs in dwarf galaxies (e.g., Rees 1988; Donley et al. 2002; 
Gezari et al. 2006; Milosavljevic et al. 2006). 
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APPENDIX 

A. MEASUREMENT UNCERTAINTIES 

We describe a suite of simulated galaxy spectra that we have built to investigate the reliability of our BH mass measurements over 
the Ha luminosity range of interest, for SDSS-quality data. We investigate our ability to recover the input FWHMh«, Lu a , and Mbh 
as the Ha flux, galaxy continuum, and S/N of the spectra change. We also justify our particular choice of detection thresholds in rms- 
normalized flux and EW(Ha). At the same time, we investigate the importance of factors such as strength of the narrow lines, shape 
of the stellar continuum, and the slope of the AGN continuum. The basic grid of simulations spans a range of 10 4 <Mbu/Mq< 10 9 , 
0.01 <Lboi/^Edd< 3, and 5 <S/N< 50. As we have noted previously, the level of galaxy continuum may play a crucial role in our 
ability to measure the broad Ha component at all. We have decided to assign the galaxy luminosity based roughly on scaling relations 
between the BH and the surrounding galaxy. While there is large scatter in such scaling relations, we attempt to span a reasonable 
range of galaxy luminosities, so that at least we are simulating realistic conditions. 

All of the spectral properties are derived from the chosen Mbh, ^boi/^Edd, and S/N. In terms of the AGN features, the continuum 
luminosity and FWHMh q are determined using the formulae presented in Greene & Ho (2005b). A small fraction (10%) of the total 
broad emission luminosity is placed in a very broad component, since such broad wings are often seen in actual spectra. The AGN 
continuum shape is given slopes of (3 = 1.5, 1.0, and 0.5, where fp oc A"* 3 . Using the Tremaine et al. (2002) Mbh-c* relation, we 
calculate the width of the narrow-line region components, assuming that the velocity dispersion of the narrow lines is equal to that 
of the bulge (e.g., Nelson & Whittle 1996; Greene & Ho 2005a), while the narrow-line flux is fixed using the empirical relation of 
Zakamska et al. (2003), which relates L[oni] to £5100- Finally, all lines are convolved to the typical velocity resolution of the SDSS, 
a ss 71 km s" 1 (e.g., Heckman et al. 2004). 

In terms of the host galaxy, we again assume the Mbh - &* relation holds, and associate a stellar velocity dispersion and bulge 
luminosity with each BH. For elliptical galaxies, this is the entire story, but at low Mbh, the hosts are typically spiral galaxies. 
Therefore, we need a way to estimate the total galaxy luminosity based on the bulge properties. The Mbh-c» relation of Tremaine 
et al. (2002) links each input BH with a <r*. Then, using the relation between cr* and maximum circular velocity v c (e.g., Ferrarese 
2002; Baes et al. 2003; Pizzella et al. 2005) to obtain v c , the Tully-Fisher relation (Tully & Fisher 1977) may be used to estimate a 
total galaxy luminosity for a spiral host. In particular, we use the best-fit Tully-Fisher relation from Masters et al. (2006) for their 
"in+" sample, and we neglect variations in the relation with Hubble type. We incorporate the scatter in both the v c -er» relation and 
the Tully-Fisher relation by perturbing each fit parameter by a log-normally distributed deviate drawn from the best-fit distribution. 
For a given BH, we must decide whether to place it in a spiral or elliptical host galaxy. There is, of course, a limit to how bright a 
spiral galaxy may be and, conversely, to the faintest elliptical galaxies. For the latter, we adopt the luminosity of M32 {Mb = -15.8 
mag; Tremaine et al. 2002), while for the former we adopt a limit of M/ > -24 mag, which represents the brightest galaxies in the 
Masters et al. sample. We note that M32 is fainter than POX 52, which is a dwarf spheroidal galaxy known to host a BH of mass 
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~ 10 5 Mq (Barth et al. 2004). This is because dwarf spheroidal galaxies are more luminous than elliptical galaxies at the same cr* 
(e.g., Geha et al. 2003). If the calculated bulge luminosity is too faint, the object is assigned only a spiral host, and conversely if the 
luminosity is bright enough, the object is assigned an elliptical host. BHs in the mass range of 10 6 5 - 10 8 M are typically assigned 
both, which is in good agreement with observations. The stellar continuum shape is provided by the eigenspectra presented in Yip et 
al. (2004). For the majority of the simulations, we use a pure absorption-line spectrum typical of an elliptical galaxy, by combining 
their first three eigenspectra with weights [1,1,0.9]. We have investigated the impact of younger stellar populations, using weights 
[0.1,0.0,-0.5] to represent a typical Sc galaxy, and [0.0,0.0,-0.5] to represent a post-starburst system. We 
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FIG. A12. — Deviation A = (out— in)/in, from the measured properties Ln a , FWHMh q , and the inferred Mbh for all simulations. We show the deviations in (left) 
logarithmic bins of Ha flux normalized by the local rms and (right) EW(Ha). The entire sample in each bin is shown as an open circle, while the sub-sample with 
the joint flux and EW cuts are shown as open squares; they only differ at and below the detection thresholds. Objects that have been rejected by our initial detection 
are not included, but those below our detection threshold are, so that we may justify our choices. The detection thresholds are indicated by dotted lines. 



find that changes in stellar populations do not significantly impact our derived masses; the PCA method of Hao et al. (2005a) is quite 
robust for our purposes. 

Each artificial spectrum is run through our entire detection algorithm. Continuum subtraction is performed using PCA, the Ha 
detection is performed, and then the line fitting is performed to see if the broad lines are detectable at this S/N. Objects that would 
be rejected by our initial Ha detection algorithm are not included in the investigation of parameters shown below, but we do not 
remove objects below our detection threshold, since it is our goal to demonstrate that we have chosen reasonable values. We quantify 
the dependence of the measurements on normalized flux and EW(Ha) by plotting as a figure of merit the fractional deviation 
AX = (X out -Xi n )/Xi n . Specifically, we calculate AL Hq , AFWHM Hct , and AM B h- Note that these deviations are plotted in linear 
(not logarithmic) space. A factor of 2 error in M Bli corresponds to M BH =1. We show in Figure A12 the dependence of AX on 
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the rms-normalized flux and EW(Ha) for each parameter. The detection threshold is noted as a dotted line in each panel. It is quite 
clear from these figures that below our chosen detection thresholds we incur completely unacceptable uncertainties in the derived 
parameters. We have run smaller sets of simulations to test the importance of both the narrow-line region strength and the underlying 
galaxy continuum shape, and we find that the deviations in L H a, FWHM Hq , and M B h are constant as a function of these parameters. 
If anything, low-mass BHs are slightly more likely to be included in our sample when the narrow-line strength is high compared to 
our fiducial choice. 



B. MAXIMUM LIKELIHOOD 

While the V /V ma x method used in §4.1 is extremely intuitive, the resulting luminosity functions may be biased due to local 
inhomogeneities in the distribution of galaxies, as well as peculiar motions within the Local Supercluster (e.g., Efstathiou et al. 
1988). In the maximum likelihood method, rather than parameterizing the AGN selection function as a maximum observable volume, 
each source is assigned a probability of detection as a function of Ha luminosity, /5,(L H a)- Then, the likelihood that we observe our 
particular sample may be expressed as 

YltPiiumLiidLi 

Our goal is to maximize this function, or equivalently minimize S = -2 In C. Rather than assume a particular form for the underlying 
luminosity function, we derive a stepwise-constant, non-parametric function. In detail, we follow the iterative procedure outlined in 
Blanton (2000; see also Koranyi & Strauss 1997; Hao et al. 2005b). 

The primary challenge is in deriving the selection function, or the probability of observing a given object as a function of lumi- 
nosity. We attempt to follow closely the technique of Hao et al. (2005b) for ease of direct comparison. Given the flux limits of the 
SDSS, each source has both a minimum and a maximum luminosity at which we may observe it. As Lhq increases, we increase the 
power-law luminosity of the AGN according to the relation of Greene & Ho (2005b), and the maximum luminosity is reached when 
the source reaches either its Eddington luminosity, or the bright limit of the SDSS. The minimum luminosity is more complicated 
to calculate. It is reached either when the source reaches a magnitude limit, or when the Ha line is no longer detectable using our 
algorithms. For objects targeted as galaxies, this is relatively straightforward to determine. At each luminosity, we calculate the 
new total magnitudes and corresponding typical S/N. The minimum luminosity is reached when the source reaches our detection 
threshold (/e a /rms < 200) or the galaxy magnitude limit. For the objects targeted as QSOs, the situation is somewhat complicated 
by the possibility that the object changes target selection from QSO to galaxy as the source fades. The AGN contribution is estimated 
from the Ha luminosity. We adopt a very ad hoc prescription, in which the object is converted from a QSO to a galaxy when the 
AGN accounts for less than 10% of the total Petrosian magnitude. In practice, the number we choose does not affect our results, 
since in the vast majority of cases the limit is set by the detection threshold. 

Within the limiting magnitudes, the AGN will have different detection thresholds depending on the contrast of the Ha luminosity 
and the S/N. As the Ha luminosity decreases for a given source, several important properties of the spectrum will change. The 
continuum level will decrease, of course, but since the BH mass is not changing, the line shape also must change. Progressively, 
as objects move to lower Eddington ratio, they grow more difficult to find, as their Ha flux is distributed more and more broadly, 
and the overall S/N of the spectrum decreases in parallel. In terms of detection efficiency, as we show in Appendix A, the galaxy 
type and luminosity only provide secondary corrections to our detection efficiency, which depends most strongly on line luminosity, 
line strength, and S/N. For that reason, we use our suite of simulated AGNs that span a wide range in M Bli , Lboi/^Edd, and S/N, as 
described in above. Our procedure is the following. We step between the minimum and maximum luminosity in steps of 0.5 dex (a 
factor of 3 in luminosity), and recalculate the total magnitude and thus typical S/N of the source. Using the measured Mbh and the 
new Ln a , we also assign an Eddington ratio to each source. All simulations with similar M sli , L\, i/L E dd, and S/N are then selected, 
and a weighted-mean detection fraction is calculated (that is, we assign more weight to the objects with the most comparable values 
in each parameter). The benefit of this approach is that the mass-dependent biases highlighted above are built into the calculation. 
We are assuming, however, that the galaxy level scales with the BH mass. If this assumption is severely violated, then the galaxy 
light may overwhelm the AGN light more rapidly than assumed in these simulations. 

Using the calculated minimum and maximum luminosities and detection fractions, we derive the maximum likelihood luminosity 
function shown in Figure 5. That we find such good agreement using two different methodologies is very encouraging. As an 
additional test, we divide our sample into redshift bins and verify that our maximum likelihood luminosity function, combined with 
the detection fractions calculated above, accurately predicts the observed distribution of L Ha in each bin (Sandage et al. 1979). 
Following Hao et al., we calculate, in each redshift bin, the predicted number of sources per L Hq bin (F ? ) as: 

F P (L)AL=V/- (L)$(L)A£ ,. (B2) 

We then compare the calculated F v with the observed L Hq distribution, as shown in Figure B13. As can be seen, the agreement is 
very satisfactory, suggesting that our luminosity function provides a good description of the data. 



BH Mass Function 
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FIG. B 13. — A comparison between the actual distribution of La a (crosses and dotted lines) and that predicted by the luminosity function and detection fractions 
(solid histograms). Each panel shows a different redshift bin. The good agreement provides an important sanity check on our luminosity function. 
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